function sys=obv1(t,x,U,D,x2)
sys=zeros(2,1);

J=2;alpha1=1;alpha2=1;lamb0=1;lamb1=0.1;lamb2=1;lamb3=0.1;m1=7/5;m2=7/5;

n1=x(1);
s=x2-n1;

sys(1) = 1/J*(U + D + alpha1*s^(2*m1-1)+lamb0*s + lamb1*sign(s));
sys(2) = lamb2*(-U-D) + alpha2*(lamb1*sign(s))^(2*m2-1)+lamb3*sign(lamb1*sign(s));
